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1.  OBJECTIVE 


1.1  In  trod uc tlon 

Metallic  elements  volatilized  during  meteor  entry  Into  the  Earth's  upper 
atmosphere  play  a  significant  role  In  the  structure  of  the  D  and  E  regions  of 
the  Ionosphere,  ~  and,  at  least  In  the  case  of  sodium,  the  visible  day  and 
nlghtglow  emissions  from  the  mesosphere  and  lower  thermosphere.^"^  Recently, 
It  was  suggested  that  sodium  and  other  meteor  metals  may  be  Important  In 
stratospheric  chemistry  by  affecting  ozone  reduction  by  the  catalytic  chlorine 
cycle . ® 

The  Influx  of  meteor  metals  Into  the  upper  atmosphere  has  been 
estimated^  to  be  3.5  x  10^  kg  yr~^,  with  a  sodium  abundance  of  2Z  leading  to  a 
calculated  sodium  flux  of  1.2  x  10^  atoms  cm~^  s~^.  Other  estimates  of  sodium 
flux  run  as  high  as  2  x  10**  cm"^  The  flux  of  otner  metallic  species 

such  as  Mg,  Ca,  Al,  SI  and  Fe  will  be  as  much  as  10  times  higher  and 
speculative  concerns  about  their  Influence  on  upper  atmospheric  homogeneous 
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and  heterogeneous  chemistry  have  been  published.  S 

Unfortunately,  all  attempts  to  model  the  role  of  volatilized  meteor 
metals  (particularly  sodium)  in  the  mesosphere  and  stratosphere^"***  ^ 
have  suffered  from  an  almost  total  lack  of  measured  rate  constants.  All  such 
models  start  with  the  oxidation  of  sodium  or  other  metallic  species  In 
reaction  with  atmospheric  0,  O2 ,  or  O3.  Until  1982,  the  only  measured 
chemical  rate  constants  available  for  any  meteor  metal  oxidation  reactions 
were  the  three  body  recombination  reactions  of  alkali  atoms  (Na,  R)  with 

O2 16-18 

Na  or  K  +  0^  +  M  -*•  NaO,  or  K0„  +  M  .  (1) 

2  2  2 

Flame  studies  by  Hynes  et  al^^  and  direct  measurements  at  724  and  844  K 
by  Husain  and  Plane showed  that  this  reaction  rate  Is  1000  times  faster  than 


1 
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previously  believed.  ~  Recent  direct  flow  tube  measurements  In  our 
laboratory  have  extended  these  results  from  700  K  down  to  300  K,  and  show  that 
the  reaction  rates  vary  Inversely  with  temperature,  and  confirm  the  larger 
values  for  the  rate  constants. 

As  an  indication  of  the  effect  that  the  availability  of  chemical  rate 
data  has  on  aeronomlc  modeling  studies,  consider  the  theoretical  calculation 
of  several  Important  sodlun  oxidation  reaction  rate  constants: 

Na  +  0^  NaO  +  0^  (2) 

NaO  +  0  ->■  Na  +  0^  (3) 


and 

NaO  +  Oj  *  Products  (4) 

performed  at  Aerodyne  Research,  Inc.^^  The  adoption  of  realistic  values  for 
just  these  three  rate  constants  has  led  to  a  major  revision  In  the 
understanding  of  the  structure  of  the  mesospheric  sodium  layer, the 
magnitude  of  the  sodium  D  line  nightglow,^ » and  the  understanding  of 
long-lived  luminous  meteor  trails.® 

From  the  recent  modeling  work^"®> and  comparison  of 
these  models  with  atmospheric  measurements,  neutral  sodium  is 

believed  to  be  transformed  via  a  series  of  chemical  reactions  Involving  NaO, 
Na02 ,  and  NaOH.  A  survey  of  the  relevant  literature^"® > provides  a 
fairly  complete  list  of  possible  neutral  sodium  reactions  of  Importance. 

These  are  listed  In  Table  1,  In  four  groups,  which  to  some  extent,  reflect  the 
sequence  of  events  from  the  Introduction  of  neutral  sodium  at  110  km  altitude 
to  Its  disappearance  near  70  km.  Also  listed  are  reaction  enthalpies,  with 
associated  uncertainties. 


NaO  +  O3  -►  NaOj  +  Oj  282  ±47 

Na  +  2O2  119  ±42 

NaO  +  H2O  -*■  NaOH  +  OH  1  ±44 

NaO  +  H2  +  NaOH  +  H  64  ±44 

NaO  +  CH4  -*•  NaOH  +  CH3  61  ±44 

NaO  +  HCl  -►  NaCl  +  OH  134  ±42 

NaO  +  0  +  Na  +  O2  225  ±42 

NaO  +  HO2  -*•  NaOH  +  O2  303  ±44 

NaO  +  H  -►  Na  +  OH  155  ±42 

NaO  +  CIO  -*•  NaOCl  +  0 


Na02  +  HCl  NaCl  +  HO2  13  ±23 

Na02  +  Cl  +  NaCl  +  O2  247  ±21 

Na02  +  CIO  ♦  NaCl  +  O3  84  ±21 

Na02  +  H2  •*>  NaOH  +  OH  103  ±24 

♦  NaO  +  H2O  103  ±^‘’ 

Na02  +  0  NaO  +  O2  110  ±47 

Na02  +  OH  NaOH  +  O2  181  ±24 

Na02  +  H  -*•  NaOH  +  0  111  ±24 

-*■  Na  +  HO2  34  ±23 

III.  Sodium  Hydroxide  Reactions 

NaOH  +  HCl  -*■  NaCl  +  H2O  133  ±13 

NaOH  +  Cl  +  NaCl  +  OH  67  ±13 

NaOH  +  CIO  NaCl  +  HO2  63  ±15 

NaOH  +  H  Na  +  H2O  154  ±13 


NaOH  +  CO2  +  NaHC03 
NaOH  +  HNO3  -»■  NaN03  +  H2O 
NaH003  +  HCl  -►  NaCl  +  Products 


IV.  Photolysis  Reactions 

NaOH  +  hv  >  Na  +  OH 

NaO  +  hv  +  Na  +  0 

Na02  +  hv  +  Na  +  O2 

NaCl  +  hv  +  Na  +  Cl 


The  main  removal  mechanisms  of  sodium  are  by  reaction  with  O2  or  O3,  the 
latter  reaction  used  by  Chapman^^  to  describe  the  Na  nlghtglow.  The  recent 
rate  measurements  of  reaction  (1)  however,  lead  us  to  believe  that  below  80  km 
most  of  the  Na  Is  converted  to  Na02. 

The  fate  of  Na02,  however,  Is  uncertain.  Recent  measurements  by  Flgger 
et  al^**  and  flame  studies  by  Hynes  et  al^^  Imply  that  the  Na-02  bond  strength 
is  much  weaker  than  previously  believed,  with  a  value  of  ca.  163  kJ  mol”^.  As 
a  result,  Na02  will  react  exothermically  with  atomic  oxygen  to  produce  NaO. 
Sodium  oxide  can  react  either  with  ozone  to  reform  Na02  or  react  with  water  to 
produce  NaOH. 

NaO^  +  0  ->•  NaO  +  0^  (5) 

NaO  +  0  -*•  NaO  +  0  (6) 

NaO  +  H^O  -*•  NaOH  +  OH  (7) 

Unfortunately,  none  of  these  rates  has  been  measured,  leading  the  various 
modelers  to  speculate  on  the  importance  of  each. As  shown  by  Sze  et 
al,  the  predominant  alkali  species  below  70  km  may  be  NaO,  Na02,  and/or 
NaOH,  depending  on  ones  choice  for  kj,  kg,  and  k7. 

O 

In  a  recent  paper  by  Murad  et  al,  it  was  proposed  that  the  reactions  of 
metal  hydroxides  and  superoxides  with  chlorine  compounds  between  40  and  70  km 
may  have  an  Impact  on  the  depletion  of  stratospheric  ozone.  In  the  case  of 
sodium,  the  exothermic  blmolecular  reactions 
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NaOH  +  HCl  ->•  NaCl  +  H^O 


(8) 


NaOH  +  CIO  -*•  NaCl  +  HO^  ,  (9) 

NaOH  +  Cl  -►  NaCl  +  OH  (10) 

NaO^  +  HCl  -y  NaCl  +  HO^  ,  (11) 

Na02  +  CIO  -y  NaCl  +  O3  ,  (12) 

NaO^  +  Cl  -y  NaCl  +  0^  (13) 


might  be  expected  to  proceed  rapidly  and  act  as  a  sink  for  Cl,  given  that  NaCl 
can  readily  polymerize  and  condense  via  heterogeneous  ni-cleation.  ®  Murad 
et  al  calculated  that  if  the  reaction  rate  of  these  sodium  species  with  HCl 
were  "10“^^  cm^  molecule”^  3“^,  they  might  be  comparable  to  the  major  Cl 
regeneration  mechanism, 


OH  +  HCl  -y  Cl  +  HjO  .  (14) 

While  previously  published  studies  have  viewed  NaCl  as  a  potential  sink  for 
stratospheric  chlorine,^”®  more  recent  analyses  by  F.S.  Rowland indicate 
that  photolysis  of  NaCl  may  In  fact  release  free  Cl.  Given  the  potentially 
large  J  values  (atmospheric  photolysis  rate)  for  this  process, reactions 
(8)-(13)  could  effectively  supplement,  rather  than  remove,  reaction  (14)  as  a 
release  mechanism  for  Cl  from  the  Inactive  HCl  stratospheric  reservoir  and 
thereby  determine  the  extent  to  which  ozone  might  be  depleted  by  chlorine 
compounds  in  the  stratosphere.  However,  to  fully  understand  the  role  of 
alkali  species  In  the  stratosphere  one  must  also  consider  the  effects  of  Na02 
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and  NaOH  photodissociation  on  those  processes.  » 


The  objective  of  this  work  Is  to  provide  direct  experimental  measurements 
of  these  key  rate  constants.  Obtaining  accurate  reaction  rates  for  the  sodium 
atmospheric  system  will  enable  more  accurate  modeling  of  the  sodium  emission 
and  density  profiles.  Once  this  is  completed,  these  results  can  be  extended 
toward  the  larger  goal  of  understanding  the  chemistry  of  the  more  abundant 
meteor  metals. 

The  work  during  this  past  year  has  concentrated  on  two  reactions,  Na  (and 
K)  +  O2  +  M,  and  NaOH  +  HCl.  Experimental  determinations  of  their  rate 
constants  are  presented  In  the  next  two  sections  of  this  report. 

1 .2  Work  Forecast 

During  the  next  year  of  this  program,  we  will  continue  our  examination  of 
the  chemistry  of  alkali  metals  In  the  mesosphere  and  upper  stratosphere.  Our 
initial  work  will  be  the  measurement  of  the  rate  constant  for  Na  with  ozone. 
This  will  also  provide  us  with  a  source  of  NaO,  with  which  we  will  measure  the 
very  Important  rate  constant  for  the  reaction  of  NaO  with  H2O.  Detection 
techniques  for  NaO  have  been  developed  during  the  first  year  of  this  program. 
We  also  plan,  later  In  the  year,  to  begin  Investigating  the  chemistry  of  NaOj. 
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2.  TEMPERATURE  DEPENDENT  TERMOLECULAR  REACTION  RATE  CONSTANTS  FOR  POTASSIUM 

AND  SODIUM  SUPEROXIDE  FORMATION 

(To  be  submitted  to  the  20*-^  Symposium  (International)  On  Combustion.) 

INTRODUCTION 

The  chemistry  of  alkali  metals  In  flames  commands  our  attention  for  both 
scholarly  and  practical  reasons.  From  a  scientific  standpoint  reactions  of 
alkali  atoms,  with  their  single  valence  electron,  form  a  natural  test  bed  to 
extend  our  theoretical  understanding  honed  on  reactions  of  atomic  hydrogen. 
Furthermore,  the  low  Ionization  potentials  of  alkali  atoms  open  the 
possibility  of  electron  transfer  reaction  mechanisms  and  allow  determination 
of  the  role  of  Ionic  potential  surfaces  In  reaction  dynamics.  Finally, 
Incredibly  sensitive  methods  of  detecting  alkali  atoms,  including  hot  wire 
surface  Ionization,  laser  Induced  resonance  fluorescence,  and  laser  Induced 
resonance  Ionization  now  allow  the  design  of  elegant  experimental  studies  for 
environments  ranging  from  molecular  beams  to  high  pressure  flames. 

From  a  practical  standpoint,  the  alkali  chemistry  of  flames  Is  focused 
primarily  on  understanding  and  utilizing  the  flame  suppression  properties  of 
alkali  additives.  Laboratory  studies in  the  1960's  and  early  1970' s 
clearly  Indicated  the  critical  nature  of  the  reaction: 

KOH  (NaOH)  +  H  ^  K  (Na)  +  H2O  (1) 

as  the  major  radical  scavenging  reaction  responsible  for  alkali  flame 
suppression.  Since  then,  detailed  computer  studies  of  alkali  suppression  of 
secondary  muzzle  flash‘d and  of  rocket  exhaust  plume  afterburning®*^  have 
confirmed  the  Importance  of  reaction  (1),  but  have  also  Indicated  that 
secondary  formation  of  alkali  superoxides: 
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K  (Na)  +  O2  +  M  >  KO2  (Na02)  +  M 


(2) 


followed  by  potential  superoxide  reactions  such  as: 

-*■  K(Na)  +  HO2 

KOj  (NaOj)  +  H  (3) 

-*>  KO(NaO)  +  OH 

KO2  (Na02)  +  H2  +  KOH  (NaOH)  +  OH  (4) 

>  KO  (NaO)  +  HO2 

KO2  (Na02)  +  OH  (5) 

-*•  KOH  (NaOH)  +  O2 

can  significantly  affect  the  degree  of  alkali  suppression  In  many 
circumstances. 

However,  for  all  the  Intrinsic  theoretical  and  experimental  Interest  In 
alkali  reactions,  our  understanding  of  the  gas  phase  oxidation  chemistry  of 
alkali  atoms  has  only  recently  progressed,  first  back  to,  and  then  beyond  the 
level  gained  In  the  1930* s  In  pioneering  sodium  diffusion  flame  studies  by 
Haber  and  Sachsse^  and  Bawn  and  Evans. ^  While  these  low  pressure  studies  of 
reaction  (2)  showed  It  proceeded  with  a  termolecular  rate  constant  In  excess 
of  10“  cm  s”  ,  >  subsequent  work  (primarily  In  the  mid  1960* s)  using 
flame  photometric  methods  by  Kaskan^®*^^  and  McEwan  and  Phillips 
Indicated  much  lower  MO2  formation  rates  (M*Na,K) as  well  as 
relatively  high  dissociation  energies  for  the  M-O2  bond  (M“Na,Ll) . 

Analyses  of  thermochemical  cycles  by  Herm  and  Herschbach^**  and 
Alexander^ ^  clearly  Indicated  that  the  values  of  M-O2  bonds  deduced  from  flame 
photometric  data  by  McEwan  and  Phillips  were  greatly  overestimated.  This  was 
recently  confirmed  by  Flgger  and  co-workers  who  analyzed  observations  of 
chemiluminescence  between  crossed  molecular  beams  of  alkali  dimers  and 
molecular  oxygen  to  yield  alkali  superoxide  bond  energies.^®  Reinterpretation 
of  the  earlier  flame  experiments^^  as  well  as  newly  performed  flame  studies 
are  also  consistent  with  lower  alkali  superoxide  bond  energies. 


The  low  sodium  and  potassium  termolecular  reaction  rates  deduced  from  the 
flame  photometry  studles^*^ » are  also  Incorrect.  Recently 
published  flash  photolysis  studies  over  limited  temperature  ranges  for  LI 
(393  and  463  K),^^  Cs  (323  K),20  Na  (724  and  844  K),21  and  K  (753  and 
853  as  well  as  flame  photometry  studies  for  all  clearly  show  fast 

termolecular  rates  for  alkali  superoxide  formation. 

In  this  paper  we  report  termolecular  rate  constants  for  superoxide 
formation  of  sodium  and  potassium  over  a  much  wider  temperature  range  than 
previously  available,  yielding  the  first  clear,  direct  measurement  of  the 
temperature  dependence  of  these  processes.  Third  body  dependencies  are 
Illustrated  by  separate  measurements  for  N2>  He  and  Ar.  These  measurements 
were  made  In  Aerodyne's  high  temperature,  fast  flow  reactor  utilizing  laser 
Induced  resonance  fluorescence  detection  of  atomic  alkali  species.  A 
discussion  of  these  results  In  terms  of  their  Importance  In  testing 
termolecular  association  reaction  theory  and  for  the  more  applied  concerns  of 
secondary  combustion  suppression  and  the  oxidation  of  meteor  deposited  metals 
In  the  upper  atmosphere  will  also  be  presented. 


EXPERIMENTAL 


The  high  temperature  fast  flow  reactor  used  In  these  experiments  has  been 
described  In  detail  prevlously^^  and  Is  shown  In  Fig.  1  In  the  configuration 
used  In  this  study.  Briefly,  a  7.26  cm  diameter,  120  cm  long  alumina  tube  Is 
heated  with  Kanthal  resistance  heaters  capable  of  maintaining  uniform 
temperatures  In  the  reaction  zone.  The  helium,  argon  or  nitrogen  carrier  gas 
enters  the  flow  tube  through  two  mulllte  multichannel  arrays  which  both 
preheat  and  lamlnarlze  the  flow  upstream  of  the  reaction  zone.  A  sufficient 
distance  (20  cm)  is  allowed  downstream  of  thes^  arrays  for  the  carrier  gas  to 
develop  a  parabolic  flow  profile  before  reaching  the  reaction  zone.  The  gas 
Is  pumped  by  a  450  liter  s“^  Roots  blower  and  mechanical  fore  pump 
combination.  Linear  flow  velocities  from  450  to  1900  cm  s~^  are  used  In  the 
present  experiments.  Carrier  gas  flow  rates  from  4U  to  400  STP  cm^  3“^  are 
measured  with  calibrated  rotameters.  Total  pressures  In  the  range  from  1  to  8 
torr  are  measured  with  a  capacitance  manometer  at  the  downstream  end  of  the 
reaction  zone.  Gas  temperatures  are  obtained  with  a  chromel-alumel 
thermocouple  which  Is  movable  throughout  the  reaction  zone.  Although  the 
apparatus  Is  capable  of  temperatures  up  to  1500  K,  experimental  considerations 
as  described  below  limited  the  upper  temperature  to  700  K  for  these  studies. 

Alkali  atoms  are  generated  by  heating  the  pure  metal  In  a  25  mm  diameter 
cylindrical  silver  plated  monel  oven  to  a  temperature  sufficient  to  obtain  a 
vapor  pressure  of  10”^  to  10“**  torr  within  the  oven.  The  oven  Is  mounted  on 
the  end  of  a  movable  13  mm  od  alumina  tube  concentric  to  the  main  flow  tube. 
The  vapor  Is  entrained  In  a  flow  of  carrier  gas  and  Introduced  Into  the  flow 
tube  In  one  of  two  configurations.  For  the  kinetic  measurements  with  O2,  the 
oven  Is  placed  at  a  fixed  position  upstream  of  the  multichannel  arrays  and  the 
alkali  vapor  Is  Introduced  into  the  main  flow  through  a  10  cm  length  of  19  mm 
o.d.  silver  tubing  which  passes  through  the  center  of  the  arrays.  In  the 
second  configuration,  the  oven  is  placed  downstream  of  the  arrays  and  is 
movable  throughout  the  reaction  zone  to  determine  the  loss  of  alkali  atoms  on 


the  reactor  walls.  The  oven  temperature  is  controlled  with  resistive  heating 
elements  independent  of  the  flow  tube  heaters.  Initial  alkali  atom 
concentrations  in  the  reaction  zone  are  maintained  at  less  than  10^*^  atoms 
cm“^  by  adjusting  the  oven  temperature  and  the  carrier  gas  flow  rate  through 
the  oven. 

The  alkali  atoms  are  detected  by  laser  Induced  fluorescence  at  the 

downstream  end  of  the  flow  tube.  A  Molectron  DL14  nitrogen  pumped  dye  laser 
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is  used  to  excite  the  4s  S1/2  -*•  5p'^P3/2  transition  at  404.4  nm  for 
potassium  or  the  3s  S1/2  3p  ^^3/2  transition  at  589.0  nm  for  sodium. 
Fluorescence  is  collected  using  a  gated  phototube  and  a  computer  controlled 
data  acquisition  system. ^3  The  signals  are  averaged  over  100  laser  pulses 
after  subtracting  for  nonfluorescent  background  contributions  and  normalizing 
for  pulse  to  pulse  fluctuations  in  laser  intensity.  Although  no  direct 
calibration  of  the  fluorescence  signal  was  attempted,,  estimates  of  sensitivity 
from  known  phototube  and  Integrator  response,  measured  laser  power,  and  atomic 
transition  probabilities  indicate  a  detection  limit  for  Na  of  10**  cm”^  and  for 
K  of  10®  cm“^  using  these  transitions. 

Molecular  oxygen  is  added  through  a  movable  loop  Injector  at  distances 
from  10  to  60  cm  upstream  of  the  detection  region.  O2  flow  rates  are  measured 
with  a  calibrated  thermal  conductivity  mass  flow  meter.  The  kinetic 
measurements  are  conducted  under  pseudo-first  order  conditions  with  2  x  10^** 
<[02]<  7  X  10^®  cm“^,  in  large  excess  compared  to  the  alkali  metal  atom 
concentration,  yet  always  <5%  of  the  total  gas  concentration.  First  order 
rate  constants  are  determined  for  a  fixed  O2  concentration  by  recording  the 
change  in  alkali  atom  concentration  as  a  function  of  O2  Injector  distance. 

Corrections  for  both  axial  and  radial  diffusion  and  wall  removal  are  made 
using  the  procedure  outlined  by  Brown.  This  method  is  based  on  a  numerical 
solution  of  the  continuity  equation  including  diffusion,  first  order  chemical 
reaction  and  wall  removal  in  a  cylindrical  cube  with  fully  developed  laminar 
flow.  A  multiplicative  correction  factor,  obtained  for  each  set  of 
experimental  conditions,  is  used  to  relate  the  observed  decay  to  the  true 
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first  order  rate  constant.  This  correction  factor  is  dependent  on  the  wall 
removal  rate  constant,  and  the  binary  diffusion  coefficient  for  the 
alkali  atom  in  the  carrier  gas.  Values  for  kw  are  determined  in  separate 
experiments  at  each  pressure  and  temperature  by  observing  the  change  in  atom 
concentration  at  the  detection  point  while  varying  the  source  oven  position. 
The  loss  of  atoms  to  the  wall  was  found  to  be  diffusion  limited  for  all  the 
conditions  used  in  this  study,  implying  a  surface  acconunoda tlon  coefficient  y 
>  0.1.  The  observed  wall  loss  rate  under  these  conditions  may  therefore  be 
used  to  determine  the  diffusion  coefficient.  This  procedure,  based  on 
Brown's^**  approach  for  obtaining  solutions  to  the  transport  equation,  is  given 
in  more  detail  elsewhere. ^5  The  resulting  values  for  D  (1  atm),  which  are 
used  in  the  diffusion  corrections  for  this  study,  are  given  in  Table  I.  The 
correction  factors  obtained  by  this  method  were  in  the  range  k^rue/kobs 
from  1.52  to  1.86,  with  moat  lying  between  1.60  and  1.65.  Under  the  flow 
conditions  of  these  experiments,  the  ratio  of  reaction  to  diffusion  times  is 
not  unlike  those  encountered  in  ion-molecule  reactions,^®  which  also  exhibit 
correction  factors  on  the  order  of  1.6. 


Table  I.  Binary  Diffusion  Coefficients  for  Na  and  K  in  He  and  N2 

2  1 

D  (cm  s~  at  1  atm) 
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Reaction  rate  constants  of  atomic  sodium  with  O2  were  measured  over  the 
temperature  range  of  320  -  698  K  for  nitrogen  as  the  third  body  collision 
partner,  and  over  the  range  of  309  -  473  K  In  helium;  for  potassium,  the 
corresponding  temperature  ranges  were  308  -  720  K  and  296  -  520  K, 
respectively.  Room  temperature  measurements  were  also  made  In  argon  for  both 
Na  and  K.  Pressures  encompassed  at  each  temperature  range  from  1  to  between  5 
and  8  torr.  Above  700  K,  no  reliable  decays  could  be  measured,  because  alkali 
atoms  which  had  accumulated  on  the  walls  at  lower  temperatures  began  to  effuse 
Into  the  flow.  Interfering  with  the  rate  measurements. 

Using  the  data  for  Na  +  O2  +  (He)  as  an  example,  we  see  In  Fig.  2  that 
the  first-order  (logarithmic)  decays  of  signal  with  reaction  time  are  linear 
for  over  a  factor  of  100.  Plots  of  a  series  of  first-order  rate  constants 
(corrected  for  wall  and  dlffuslonal  effects)  versus  [O2]  are  Illustrated  In 
Fig.  3.  Second-order  rate  constants  derived  from  these  are  then  plotted  vs. 
total  number  density  (Fig.  4).  The  third-order  rate  constants  are  then 
determined  as  the  slopes  of  these  plots.  In  most  cases,  these  lines  Intersect 
the  origin  to  within  the  statistical  uncertainties  of  the  fit.  For  Na  +  O2  + 
(N2),  a  finite  Intercept  Is  observed  at  all  temperatures,  although  It  never 
exceeds  4  x  10“^^  cm^  s“^.  This  may  be  attributed  to  the  fact  that  the  oxygen 
concentration  reached  10-15X  of  the  total  number  density  In  the  low  pressure 
runs  for  this  system.  If  kIII(M  »  O2)  >  k^^^(M  =  N2),  then  what  Is 
reported  as  klI(N2)  would  be  overestimated,  resulting  In  a  small  finite 
Intercept  and  a  value  for  k^^^(Na  +  O2  +  N2)  that  Is  ~  15%  low.  An 
alternating  explanation  sometimes  given  for  such  behavior  Is  that  the  wall 
acts  as  a  stabilizing  third  body  and  adds  a  pressure  Independent  component  to 
kH.  However,  this  explanation  would  appear  unlikely  In  these  experiments 
given  the  near  unit  efficiency  for  wall  rerr’oval  of  all  alkali  species. 

The  three-body  association  rate  constants  as  a  function  of  temperature 
for  various  third  body  collision  partners  are  shown  for  sodium  In  Fig.  5  and 


n 


T.jb  R  A  2.  k  R  ^  »  . 


for  potassium  In  Fig.  6.  The  data  can  be  described  by  in  expression  of  the 
form  klll(T)  =  AT“^.  The  results  of  nonlinear  least-square  fits  to  this 
function  are  found  in  Tables  II  and  III.  For  sodium  in  helium  and  nitrogen, 
and  for  potassium  in  helium,  a  simple  T“^  representation  could  also  be  used. 
However,  the  results  for  potassium  in  nitrogen  show  a  weaker  temperature 
dependence,  best  represented  by  ± 

The  experimental  values  for  these  rate  constants  are  listed  in  Tables  II 
and  III.  The  expressed  total  experimental  uncertainties,  including  allowance 
for  systematic  errors,  can  be  estimated  as  the  square  root  of  the  sum  of  the 
squared  individual  uncertainties  due  to  a)  flow  velocity,  temperature, 
pressure,  and  concentration  calibration  factors,  5%,  b)  the  precision  in 
determining  the  coefficient  A  in  the  expression  for  klll(T),  17%,  and  c)  the 
uncertainty  in  the  correction  of  observed  first-order  decays  for  wall  removal 
and  diffusional  effects,  10%. 

Tables  II  and  III  also  contain  a  comparison  with  other  determinations  of 
k2.  Husain  and  Plane^^»^^  used  a  flash  photolysis  system  to  generate  alkali 
atoms  from  KI  and  Nal  and  followed  the  decay  of  R  and  Na  by  atomic  resonance 
absorption.  They  observed  no  temperature  dependence  over  their  limited 
temperature  range  (724  R  and  844  K)  for  the  Na  +  O2  +  M  reaction.  A 
substantially  larger  variation  with  temperature  was  obtained  for  their  R  +  Oj 
+  M  study  although  they  arbitrarily  fit  their  data  to  a  T~^  dependence.  The 
agreement  between  their  values  and  the  slightly  extrapolated  temperature 
dependent  rate  constants  for  both  reactions  from  this  work  is  satisfactory, 
especially  considering  the  very  different  methods  employed.  This  comparison 
is  shown  in  Figs.  5  and  6. 
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Hynes  et  al.  have  studied  the  Na  +  O2  +  M  reaction  in  experiments  using 
laser  Induced  fluorescence  to  measure  Na  and  OH  profiles  under  fuel  lean 
conditions  in  H2/O2/N2  flames.  A  detailed  kinetic  model  assuming  reasonable 
values  for  rate  constants  for  the  reactions  interconnecting  the  species  Na, 
Na02 ,  NaO  and  NaOH  is  used  to  determine  the  best  fit  to  the  observations. 

Their  "best  fit”  value  of  k(Na  +  O2  +  M),  for  which  they  state  an  accuracy  of 
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a  factor  of  two,  Is  reported  assuming  a  T~^  temperature  dependence.  As  shown 
In  Fig.  5,  their  value  Is  lower  by  about  a  factor  of  4  compared  to  the 
extrapolated  value  from  this  work.  Hynes  et  al.^®  also  explain  the 
discrepancy  between  their  result  and  the  earlier  flame  studies, which 
obtained  a  lower  value  for  k(Na  +  O2  +  M)  by  nearly  three  orders  of  magnitude, 
as  being  due  to  the  neglect  of  the  effects  of  secondary  reactions  In  the 
Interpretation  of  the  earlier  studies.  A  similar  explanation  has  been 
proposed  by  Jensen. 

The  much  earlier  diffusion  flame  studies  *  correctly  Identified  the  Na 
+  O2  +  M  reaction  to  be  fast  for  a  three  body  process.  The  later  and  more 
extensive  study  by  Bawn  and  Evans®  obtained  a  value  at  533  K  which  Is  higher 
by  a  factor  of  3  than  this  study.  They  also  observed  a  fall-off  behavior  for 
their  second  order  rate  constant  with  pressure  above  about  10  torr  and  present 
their  results  in  terms  of  an  energy  transfer  mechanism  using  a  Llndemann-type 
expression  to  Interpret  their  data.  This  fall-off  behavior  was  not  observed 
by  Husain  and  Plane^^  who  show  linear  plots  of  k2  versus  pressure  from  25  to 
150  torr.  Their  effective  bimolecular  rate  at  150  torr  and  724  K  for  the  Na  + 
O2  +  (N2)  reaction  of  2.2  x  10“ cm®  s“^  Is  nearly  equal  to  the  high  pressure 
limit  inferred  by  Bawn  and  Evans  of  3  x  10"  cm  s“  .  It  thus  appears  that 
the  pressure  dependence  observed  by  Bawn  and  Evans  may  have  been  an  artifact 
of  their  experimental  diffusion  flame  technique.  No  curvature  was  observed  In 
our  lower  pressure  P  £  7  torr  conditions  of  this  study  as  shown  In  Fig.  4, 
and,  under  these  conditions,  the  recombination  reaction  appears  to  be  In  Its 
low  pressure,  third  order,  limit. 

In  this  limit,  the  recombination  rate  constant  can  best  be  calculated 

using  unlmolecular  rate  theory^^.  In  which  second-order  rate  constants  for 

unlmolecular  decomposition  are  related  through  the  equilibrium  constant  to  the 

third  order  (low  pressure)  association  rate  constants.  The  unlmolecular  rate 

constants  are  obtained  using  simplified  equations  of  Troe  based  on  RRKM 

theory.^® The  recombination  rate  constant  can  be  expressed  as  the 

sc 

product  of  a  strong  collision  rate  constant  k^^^  ^  and  a  weak  collision 
deactivation  efficiency  term  0^., 
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(6) 


k  =•  k  p 

rec  rec,o 

In  effect,  P^.  Is  a  terra  expressing  the  efficiency  by  which  (Na02)*  is 
stabilized  upon  collision.  In  the  absence  of  detailed  state- to- state  energy 
transfer  rates,  p^  is  an  adjustable  parameter  which  depends  on  the  average 
change  in  internal  energy  of  the  transition  state  complex  per  collision,  <AE>, 
so  that 


(7) 


Fg  is  a  correction  term  for  the  energy  dependence  of  the  density  of 
states  in  the  transition  complex.  Using  structural  data  from  Ref.  15,  and 
deriving  »  -55.7  kJ  mole“^  from  a  value  of  163  kJ  mole“^  for  the 
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Na-02  bond  energy,  Patrick  and  Golden  have  calculated  krec  Na  +  O2  + 

rt  s  c 

(He)  at  300  and  700  The  computations  result  in  k  (300  K)  *• 

L  6C  ^  O 

7.8  X  10“^®  cm®  molecule”  s“^  for  M  »  He.  Thus,  a  value  of  3^  *  0.18  is 
required  to  reproduce  the  experimental  result.  The  average  energy  transferred 
per  collision  is  0.8  kJ  mole”^ .  At  700  K,  k**^  »  4.2  x  10"^®  cm^  molecule”^ 

s~^  ,  suggesting  a  value  of  Pc  ”  0.09  if  <AE>  Is  assumed  to  be  independent  of 
temperature.  Extrapolation  of  the  experimental  value  for  helium  to  700  K 
results  in  an  observed  Pc  =•  0.14,  in  good  agreement,  considering  the 
simplified  assumptions  of  this  approach. 

This  unlmolecular  approach  has  been  shown  to  be  fairly  reliable  in 
calculating  termolecular  association  rate  constants  for  a  variety  of 
species.  >  The  values  of  Pc  M  *  N2  and  Ar  tend  to  lie  in  the  range 

of  0.1  -  0.5.  The  calculations  for  Na  +  O2  in  helium  agree  quite  well  with 
these  numbers,  given  that  one  might  expect  helium  to  be  somewhat  less 
efficient  than  N2  in  quenching  the  excited  intermediate.  For  the  similar 
reaction  H  +  O2  +  (N2),  the  calculations  also  perform  well,  giving  Pc  * 
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0.08,  even  though  the  rate  constant  for  this  reaction,  5.9  x  10“^^  cm^  s“^  at 
300  Is  approximately  30  times  slower  than  the  corresponding  sodium 

reaction. 

In  extrapolating  our  measured  rate  constants  to  higher  pressures,  it  is 
Important  to  understand  the  fall-off  behavior  as  the  reaction  mechanism  goes 
from  Its  low  pressure  third  order  limit  to  Its  high  pressure  second  order 
limit.  A  method  for  calculating  the  fall-off  parameter  based  on  RRKM  theory 
has  been  developed  by  Luther  and  Troe.^^  Using  their  approach,  the 
blmolecular  rate  constant  In  the  transition  region  may  be  calculated  as^^»^^ 

,11  I*-  "ot"] 

"  ■  k  ^  k  [Hi 

oo  o'- 

where  k®  and  are  the  high  pressure  and  low  pressure  11ml t’ng  rate 
constants,  respectively.  The  broadening  parameter,  F,  is  dependent  on  the 

molecular  structure  of  the  adduct  and  typically  Is  on  Ci.e  order  of  0.6  for 
2  7 

small  molecules.  The  Patrick  and  Golden  calculation  for  Na02*  gives  F  =■ 
0.5.^° 

In  order  to  assess  the  fall-off  behavior,  the  association  rate  k® 
may  be  estimated  from  the  dynamics  of  the  Na  +  O2  encounter.  Alkali  atom 
reactions  are  often  described  In  terms  of  their  Ionic  character  which  leads  to 
very  fast  reaction  rates  via  an  electron  jump  mechanism. Ionic  forces  In  the 
alkali-oxygen  associations,  however,  are  not  as  dominant  as  In  other  alkali 
atom  reaction  mechanisms.  Although  the  approaching  species  do  adlabatlcally 
transfer  from  the  Incident  covalent  potential  surface  to  an  Ionic  surface 
(electron  jump  mechanism),  the  crossing  distance  (r^)  is  relatively  small  as 
compared  with  alkali-halogen  systems  and  the  overall  association  rate  Is 
dominated  by  the  dispersion  Interaction  on  Che  covalent  potential  surface. 

This  Is  Illustrated  by  comparing  the  rate  constant  based  on  the  product  of  Che 
mean  collision  velocity  and  the  electron  jump  cross  section  (modified  to 
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account  for  the  dispersion  and  angular  momentum  terms  in  the  covalent 
potential  function) ,  to  a  rate  constant  computed  solely  from  a  close 
collision  model, in  which  all  collisions  surmounting  the  angular  momentum 
barrier  of  a  Cr~^  potential  contribute  to  the  rate.  For  Na  +  O2  at  300  K,  an 
electron  jump  rate  constant  (r^  =  2.6  A)  Is  1.5  x  10“^*^  cm^  s“^,  while  the 
close  collision  constant  Is  5.9  x  10“^®  cm^  s“^.  Thus,  It  appears  that  the 
rate  of  intermediate  complex  formation  for  reaction  (2)  Is  governed  by  long 
range  dispersion  forces,  as  found  In  most  covalent  systems.  Once  formed, 
however,  these  intermediates  are  certainly  (alkali)'*‘02"  Ion  pairs.  In 
either  case,  the  association  rate  constant  is  expected  to  be  fast  with  a  value 
for  koD  >  10"^^  cm^  s”^. 

Using  a  value  for  k^,  =  2  x  10“^^  cm^  s~^,  Eq.  (8)  may  be  used  to 
estimate  the  ratio  of  the  observed  rate  constant,  kQt)s>  bo  the  true  third 
order  low  pressure  limit,  k^.  At  the  maximum  pressure  In  odr  studies,  7 
torr,  one  obtains  kQ^g/ko  *  0.91  which  Indicates  that  our  experiments  are 
essentially  in  the  low  pressure  limit.  For  the  conditions  of  Husain  and 
Plane^^  (150  torr,  724  K) ,  this  same  calculation  would  predict  k^^jg/ko  = 

0.87.  For  the  combustion  conditions  with  P  ■  760  torr  T  ■  2000  K  applicable 
to  the  Hynes  et  al.  study^®,  one  obtains  kQjjg/ko  “  0.89.  Thus,  the  large 
value  for  k®  would  Indicate  that  both  of  these  studies  are  also  essentially 
In  the  low  pressure  limit  with  regard  to  the  recombination  reaction. 

Implications  for  Secondary  Combustion  Suppression 

Solid  propellants  for  both  tactical  rockets  and  artillery  shells  are 
frequently  loaded  with  potassium  salts  (or  less  frequently  sodium  salts)  to 
suppress  the  secondary  combustion  of  the  high  levels  of  H2  and  CO  found  in  the 
exhaust  gases  of  these  systems.  This  secondary  combustion  is  termed 
afterburning  for  rocket  systems  and  often  leads  to  bright  infrared®  and 
vlsible/ultravlolet^®  signatures  which  are  undesirable  In  weapons  systems. 
Secondary  combustion  in  gun  systems  is  termed  secondary  muzzle  flash,  and  also 
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produces  undesirable  weapons  signatures  as  well  as  acoustic  overpressures 
which  can  damage  gun  crew  hearing. 

Propellant  designers  are  interested  In  predicting  the  minimum  levels  of 
potassium  salt  loadings  necessary  to  produce  enough  KOH  in  the  primary 
propellant  combustion  to  reliably  suppress  secondary  combustion,  through 
reaction  (1): 


KOH  (NaOH)  +  H  t  K  (Na)  +  H2O  (1) 

when  fuel  rich  exhaust  gases  mix  with  the  atmosphere.  As  noted  in  the 
introduction,  atomic  potassium  or  sodium,  formed  either  by  reaction  (1)  or  In 
the  primary  combustion  process  will  rapidly  form  alkali  superoxide  molecules 
via  reaction  (2)  at  the  fast  rates  measured  in  this  work.  Computer  models  of 
this  process  for  both  rocket  exhaust  afterburning^"^  and  secondary  muzzle 
flash  from  guns‘*~^  indicate  that  this  process  will  be  particularly  effective 
in  the  cooler,  oxygen  rich  regions  created  by  substant  j1  atmospheric  dilution 
of  exhaust  products.  These  models  also  indicate  that,  if  the  alkali 
superoxide  formed  in  these  regions  can  react  with  residual  flame  radicals  or 
molecular  hydrogen  via  reactions  (3-5),  there  will  be  a  significant 
modification  (a  reduction  for  most  reasonable  rate  assumptions)  of  the  flame 
suppression  characteristics  of  KOH.^"^  The  large  rate  constants  and  moderate 
temperature  dependence  measured  in  this  work  for  superoxide  formation 
Increase  the  potential  role  of  alkali  superoxide  reactions  in  alkali  flame 
suppression.  Efforts  to  provide  the  first  experimental  determination  of  rate 
constants  for  alkali  superoxide  reactions  (3-5)  are  now  underway  in  our 
labora  tory . 

Oxidation  of  Meteor  Metals  in  the  Upper  Atmosphere 

Atomic  sodium  ablated  from  meteors  in  the  90-110  km  altitude  region  is 
estimated  to  be  equivalent  to  a  time  averaged  flux  of  over  10  atoms  cm"'^ 
sec"^.^^"^^  Lesser  amounts  of  LI  and  K  are  also  known  to  be  deposited  In  the 


same  manner.  Recent  Insight  Into  the  chemical  dynamics  of  the  upper 
atmospheric  sodium/oxygen  system^^  helped  stimulate  a  revival  of  Interest  in 
the  Chapman  mechanism  for  the  sodium  alrglow.***^  The  formation  of  sodium 
superoxide  through  reaction  (2),  is  now  recognized  as  the  lower  altitude 
cut-off  sink  for  free  sodium  in  the  70-80  km  region.  The  large  rate  constants 
for  reaction  (2)  reported  in  this  work  and  reference  21,  coupled  with  accurate 
lidar  profiles  of  mesospheric  atomic  sodium,  pose  a  stringent  test  for  models 
of  diurnal  and  seasonal  mesospheric  variations, and  have  caused  the 
revision  of  a  number  of  previous  studies. 

Below  90  km  a  series  of  chemical  and  photochemical  reactions  are  thought 
to  partition  atmospheric  sodium  among  Na02,  NaO,  NaOH,  and  NaCl. 

During  daylight^  NaOH****  and  NaCl**^  are  both  predicted  to  undergo  relatively 
rapid  photolysis,  releasing  atomic  sodium  in  the  upper  stratospheric  and  lower 
mesosphere.  The  large  rate  constant  measured  for  reaction  (2)  in  this  work 
ensures  that  rapid  Na02  formation  is  the  fate  of  the  overwhelming  proportion 
of  the  atomic  sodium  thus  released,  and  that  nearly  all  the  sodium  present  in 
gas  phase  molecular  form  is  routinely  cycled  through  the  Na02  molecule. 

Efforts  are  now  underway  in  our  laboratory  to  measure  reaction  rates  for  Na02 
with  labile  upper  atmospheric  species  as  well  as  the  Na02  photolysis  rate. 
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FIGURE  CAPTIONS 


Figure  1.  Schematic  view  of  the  high  temperature  flow  reactor. 

Figure  2.  Pseudo-first  order  decay  of  Na  fluorescence  signal  vs.  reaction 
distance  for  M  »  He  at  309  K  and  5.05  torr.  [O2]  “  1*09  (a), 

3.36  (b),  7.29  (c) ,  14.3  (d),  26.6  (e);  units  of  10^**  cm-\ 

Figure  3.  Plots  of  first-order  Na  decay  rates  vs.  [02]  for  M  ■  He  and  T  = 
309-335  K.  Total  pressure  (torr)  -  1.05  (o),  1.51  (o),  2.52 
(+) ,  3.04  (A)t  5.05  (7).  Each  point  has  been  corrected  for  the 
effects  of  diffusion  and  wall  loss. 

Figure  4.  Plot  of  second-order  Na  rate  constant  vs.  total  number  density  for 
M  -  He  at  <T>  -  322  K  (o)  and  472  K  (o). 

Figure  5.  Plot  of  In  k(T)  vs.  In  T  for  the  reaction  Na  +  O2  +  M,  with  M  =  N2 
(circles),  He  (triangles),  and  Ar  (square)  for  this  work.  The 
solid  lines  are  least-square  fits  to  these  data.  The  +  and  x 
symbols  represent  data  for  M  -  Nj  and  He,  respectively  from 
reference  21,  and  the  dashed  line  (M  -  flame  gases)  Is  from 
reference  18. 

Figure  6.  Plot  of  In  k(T)  vs.  In  T  for  the  reaction  K  +  O2  +  M,  with  M  »  N2 
(circles).  He  (triangles),  and  Ar  (square).  The  solid  lines  are 
least-square  fits  to  these  data.  The  -f  and  x  symbols  represent 
data  for  M  =  N2  and  He,  respectively  from  reference  22. 
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FIGURE 


3.  THE  GAS  PHASE  REACTION  RATE  OF  SODIUM  HYDROXIDE  WITH  HYDROCHLORIC  ACID 


(To  be  submitted  to  the  Journal  of  Physical  Chemistry.) 

ABSTRACT 

The  reactions  of  metallic  species  Introduced  Into  the  atmosphere  by 
meteor  ablation  may  play  a  significant  role  in  stratospheric  chemistry.  In 
particular,  it  has  been  suggested  that  the  reaction  of  NaOH  with  HCl  might 
affect  the  concentration  of  odd  chlorine,  thus  having  an  impact  on  the  ozone 
balance.  This  paper  describes  the  first  measurement  of  this  reaction  rate 
constant.  At  308  K,  we  find  that  k  =  (2.8  ±0.9)  x  10“^°  cm^  molecule"^  s“^. 

As  a  result  of  the  methods  developed  to  perform  this  measurement,  we  have  also 
determined  estimates  of  the  following  room  temperature  rate  constants  In  units 
of  cm^  molecule"^  s“^;  k(NaO  +  HCl  NaCl  +  OH)  =  2.8  x  10"^°,  k(NaOH  +  H 
Na  +•  H2O)  >  4  X  10"^^,  k(NaCl  +  H  >  Na  +  HCl)  «  5  x  10“^^^  and  k(Na  + 

H2O2)  =  6.9  X  10“^'',  where  approximately  0.6  of  the  reactions  produce  NaOH  + 
OH,  with  the  remainder  forming  NaO  +  H2O. 

INTRODUCTION 

Metallic  elements  volatilized  during  meteor  entry  into  the  Earth's  upper 
atmosphere  play  a  significant  role  in  the  structure  of  the  D  and  E  regions  of 
the  Ionosphere, and,  at  least  In  the  case  of  sodium,  the  visible  day  and 
nightglow  emissions  from  the  mesosphere  and  lower  thermosphere.^"^  Recently, 
it  was  suggested  that  sodium  and  other  meteor  metals  may  be  important  in 
stratospheric  chemistry  by  affecting  ozone  reduction  by  the  catalytic  chlorine 
cycle. 

The  Influx  of  meteor  metals  into  the  upper  atmosphere  has  been  estimated^ 
to  be  3.5  X  10^  kg  yr"^,  with  a  sodium  abundance  of  2%  leading  to  a  calculated 
sodium  flux  of  1.2  X  10**  atoms  cm“^  s“^.  Other  estimates  of  sodium  flux  run 


as  high  as  2  x  10*^  cm"^  The  flux  of  other  metallic  species  such  as 

Mg,  Ca,  Al,  SI  and  Fe  will  be  as  much  as  10  times  higher  and  speculative 
concerns  about  their  Influence  on  upper  atmospheric  homogeneous  and 
heterogeneous  chemistry  have  been  published.  ^ 

Unfortunately,  all  attempts  to  model  the  role  of  volatilized  meteor 
metals  (particularly  sodium)  In  the  mesosphere  and  stratosphere^""^’^ 
have  suffered  from  an  almost  total  lack  of  measured  rate  constants.  All  such 
models  start  with  the  oxidation  of  sodium  or  other  metallic  species  In 
reaction  with  atmospheric  0,  O2 ,  or  O3 .  However,  the  only  measured  chemical 
rate  constants  currently  available  for  any  alkali  meteor  metal  oxidation 
reactions  are  those  for  the  three  body  recombination  reactions  of  alkali  atoms 
(Na,  K)  with 

Na  or  K  +  0^  +  M  NaO^  or  KO^  +  M  .  (1) 

In  fact,  until  two  years  ago,  even  this  reaction  was  thought  to  be  ~1000 
19  20 

times  slower  “  than  the  value  found  In  the  more  recent  measurements. 
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From  the  recent  modeling  work  comparison  of 

these  models  with  atmospheric  measurements ,  ^  ^  neutral  sodium  Is 

believed  to  be  transformed  via  a  series  of  chemical  reactions  Involving  NaO 
and  Na02  Intermediates  to  NaOH.  A  survey  of  the  relevant 
literature^"®  >  provides  a  fairly  complete  list  of  possible  neutral 

sodium  reactions  of  Importance,  and  a  simplified  schematic  diagram  of  these 
reactions  Is  shown  In  Figure  1. 

Neutral  sodium  appears  In  a  range  between  110  and  70  km  altitude.  As 
shown  In  Figure  1,  the  main  removal  mechanisms  are  by  reaction  with  O2  or  O3, 
the  latter  used  by  Chapman^^  to  explain  the  Na  nlghtglow.  We  have  recently 
completed  temperature  dependent  measurements  of  the  reaction  rates  with  O2  In 
our  laboratory  which  show  that  this  reaction  proceeds  with  a  fast  three  body 
rate  constant  of  1.9  x  10"®°  cm®  s"®  with  N2  as  the  third  body  at  300  K. 

This,  leads  us  to  believe  that  most  of  the  Na  In  the  lower  mesosphere  Is 
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converted  to  Na02>  since  the  three  body  recombination  rate  with  O2  exceeds 
even  a  gas  kinetic  two  body  rate  with  O3  below  80  km.  Since  none  of  the 
remaining  rates  have  ever  been  measured,  the  rest  of  the  mechanism  (as 
proposed  by  various  models)  is  speculative.  However,  the  dominant  sodium 
species  which  leave  the  mesosphere  and  enter  the  stratosphere  are  probably 
NaOH  and/or  NaOj. 

The  fate  of  NaOH  is  very  uncertain.  Indeed,  only  in  the  past  few  years 
has  NaOH  been  recognized  as  a  major  reaction  product.  Ferguson^  suggested 
that  it  forms  NaOH  cluster  ions  of  the  form  H+(NaOH)n  '^^ich  either 

may  be  rained  out  from  the  troposphere  or  removed  in  the  stratospheric  aerosol 
layer.  In  a  recent  paper  by  Murad  et  al,®  it  was  proposed  that  the  reactions 
of  metal  hydroxides  (and  superoxides)  with  chlorine  compounds  between  40  and 
70  km  may  have  an  impact  on  the  depletion  of  stratospheric  ozone.  In  the  case 
of  sodium,  the  exothermic  blmolecular  reactions 

NaOH  +  HCl  -*■  NaCl  +  H_0  AH°  =  -133  ±13  kJ  mol"^  ,  (2) 

NaOH  +  CIO  -<•  NaCl  +  HO^  AH°  =  -63  ±15  kJ  mol”^  ,  and  (3) 

NaOH  +  Cl  NaCl  +  OH  AH®  =  -67  ±13  kJ  mol"^  (4) 

r 

might  be  expected  to  proceed  rapidly  and  act  as  a  sink  for  Cl,  given  that  NaCl 
can  readily  polymerize  and  condense  via  heterogeneous  nucleation.®  Murad 
et  al  calculated  that  if  k2  were  -lO-^’  cm^  molecule” *  s"'',  then  the 
reaction  of  NaOH  with  HCl  would  be  comparable  to  the  major  recognized  Cl 
regeneration  mechanism, 


3-^ 


OH  +  HCl  Cl  +  H  0 
2 


(5) 


While  previously  published  studies  have  viewed  NaCl  as  a  potential  sink  for 
stratospheric  chlorine,^”®  more  recent  analyses  by  F.  Rowland 2**  Indicate  that 
photolysis  of  NaCl  may  In  fact  release  free  Cl.  Given  the  potentially  large  J 
values  (photolysis  rates)  for  this  process,^**  reactions  (2)  -  (4)  could 
effectively  supplement,  rather  than  remove,  reaction  (5)  as  a  release 
mechanism  for  Cl  from  the  Inactive  HCl  stratospheric  reservoir  and  thereby 
determine  the  extent  to  which  ozone  might  be  depleted  by  chlorine  compounds  In 
the  stratosphere.  However,  to  understand  fully  the  role  of  alkali  species  In 
the  atmosphere,  one  must  also  consider  the  effects  of  Na02  and  NaOH 
photodlssoclatlon  on  these  processes.^ 

The  purpose  of  this  study  is  to  provide  a  direct  experimental  rate 
measurement  of  the  reaction  NaOH  +  HCl  -►  NaCl  +  H2O,  which  serves  as  a 
starting  point  for  understanding  the  stratospheric  role  of  sodlun  and  other 
meteor  metals.  Only  by  obtaining  accurate,  directly  measured  rate  constants 
can  we  hope  to  understand  atmospheric  metal  chemistry,  and  In  particular,  how 
this  chemistry  affects  the  ozone  balance. 


EXPERIMENTAL 

Very  little  gas  phase  research  has  been  done  on  alkali  hydroxides  because 
they  are  difficult  to  vaporize, are  extremely  corrosive, and  readily 
dimerize  In  the  gas  phase.  In  addition,  there  are  no  established  detection 
techniques  sensitive  enough  to  make  kinetic  analyses  of  alkali  hydroxide 
reactions.  With  this  In  mind  we  have  taken  an  indirect  approach  for  producing 
and  detecting  NaOH,  l.e.,  chemical  production  of  NaOH  by  the  reaction  of 
atomic  sodium  with  hydrogen  peroxide,  and  detection  of  NaOH  by  chemical 
conversion  back  to  atomic  sodium,  which  Is  then  observed  using  laser  Induced 
fluorescence. 

The  measurements  were  performed  In  the  Aerodyne  high  temperature  fast 
flow  reactor,  which  Is  fully  described  in  Reference  29  and  whose  relevant 
features  are  illustrated  In  Figure  2.  Briefly,  a  7.26  cm  diameter,  120  cm 
long  alumina  tube  Is  used,  fitted  with  Kanthal  heater  elements  which  can 


radlatlvely  heat  the  tube  over  the  temperature  range  294  K  -  1500  K.  Four 
perpendicular  alumina  side  arms  at  the  tube  exit  permit  detection  of  flow 
species  by  a  variety  of  optical  techniques,  Including  laser-  or  resonance 
lamp-excited  fluorescence,  Infrared  absorption  (either  broadband  or  high 
resolution),  and  observation  of  chemiluminescence.  Gas  temperatures  are 
obtained  with  chromel-alumel  or  shielded  W-5%Re/W-26%Re  thermocouples. 
Extensive  calibrations  Indicate  the  temperature  can  be  measured  with  an 
accuracy  of  ±10  K  over  the  entire  operating  range  of  the  reactor  for  flow 
Reynolds  numbers  below  500;  for  the  current  experiment  Is  typically  below 
50.  The  gas  Is  pumped  by  a  Kinney  KMBD  1602  mechanical  pump  and  Roots  blower 
with  an  effective  pumping  speed  of  450  liters  s“^.  The  helium  carrier  gas  Is 
added  at  the  entrance  of  the  flow  tube  through  mulllte  multichannel  arrays 
which  lamlnarlze  the  flow.  These  are  ~2.5  cm  upstream  of  the  outlet  of  the 
reactant  Inlet  tube  and  78  cm  from  the  detection  region.  This  allows 
sufficient  distance  <18  cm)  for  the  helium  flow  to  mix  with  the  reactants  and 
develop  a  parabolic  velocity  profile  before  reaching  the  reaction  zone.  Gas 
volumetric  flow  rates  are  determined  using  calibrated  thermal  conductivity 
type  mass  flow  meters.  Flow  speeds  can  be  varied  from  4  to  100  m  s~^.  A 
calibrated  MKS  Baratron  Model  310-BSH10  capacitance  manometer  (0.8%  accuracy) 
Is  used  to  measure  pressure. 

Alkali  atoms  are  generated  by  heating  the  sample  In  a  2.5  cm  dia 
cylindrical  monel  oven  to  a  temperature  commensurate  with  attaining  a  vapor 
pressure  of  that  species  of  10”^  to  10~"^  torr  within  the  oven.  The  oven  Is 
silver  plated  to  resist  alkali  corrosion.  The  vapor  Is  entrat.ied  In  a  flow 
of  Inert  carrier  gas  and  Introduced  Into  the  flow  tube  through  a  10  cm  section 
of  19  nm  o.d.  silver  tubing.  The  sodium  vapor  is  further  diluted  by  the 
carrier  gas  In  the  main  flow  tube  so  that  the  sodium  concentration  within  the 
reaction  zone  Is  always  less  than  10^*^  cm“^.  For  wall  removal  rate 
measurements  the  entire  oven  assembly  may  be  placed  downstream  of  the  mulllte 
arrays  as  a  movable  source.  Since  the  oven  is  heated.  It  warms  the  main 
carrier  gas  flow  slightly.  Axial  temperature  surveys  In  the  reaction  region 
show  that  the  final  flew  temperature  profile  Is  uniform  at  a  value  of  308  K. 


Sodium  hydroxide  is  produced  via  the  reaction 


Na  +  H2O2  ->■  NaOH  +  OH 


AH°  =•  -  130  ±13  kJ  mol"^ 


(6) 


The  H2O2  is  added  upstream  of  the  mullite  arrays  with  the  main  carrier  gas. 

An  all  glass  and  teflon  inlet  system  is  used  to  prevent  decomposition  of  H2O2 
between  its  reservoir  and  the  flow  tube.  Flow  rates  of  H2O2  are  determined  by 
entraining  the  vapor  in  a  measured  helium  flow  at  a  known  total  reservoir 
pressure,  where 


2  2 

Flow  =  (-p - )  X  (flow  helium  through  system)  (7) 

^  ^total 

The  partial  pressure  of  H2O2  at  25‘’C  is  2.0  torr. 

Detection  of  NaOH  is  accomplished  by  converting  it  back  to  atomic  sodium 
in  the  detection  region,  where  the  sodium  is  measured  by  laser  induced 
fluorescence  (LIF).  This  conversion  is  accomplished  by  injecting  an  excess  of 
atomic  hydrogen  into  the  flow  2  cm  upstream  of  the  LIF  detector. 

NaOH  +  H  ->•  Na  +  H2O  AH°  -  -154  ±13  kJ  mol~^  (8) 

Na  +  hv  -*■  Na*  -*•  Na  +  hv  LIF  Detection  (9) 

The  hydrogen  atoms,  produced  by  microwave  discharge  of  pure  molecular 

hydrogen,  are  Introduced  through  a  teflon-lined  6  mm  dia.  tube.  With  the  mean 

flow  velocity  of  1000  cm  s“^ ,  the  reaction  time  of  hydrogen  in  the  detector 

( tfj)  is  2  ms.  The  atomic  hydrogen  concentration  is  estimated  to  be  2  x  10^“* 

cm“^ ,  by  measuring  the  H2  flow  rate  and  assuming  that  10%  of  the  H2  passing 

through  the  discharge  dissociates,^^  with  no  recombination^^  or  loss  on  the 
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inlet  tube  walls. 


A  Molectron  DL14  nitrogen  pumped  dye  laser  Is  used  for  laser  Induced 
fluorescence  detection  of  Na.  The  laser  Induced  fluorescence  data  acquisition 
system,  and  manipulation  of  data  have  been  detailed  elsewhere. However,  it 
should  be  noted  that  fluorescence  is  usually  averaged  over  100  laser  pulses, 
accounting  for  nonfluorescent  background  signals  and  for  pulse-to-  pulse 
fluctuations  In  laser  Intensity.  The  combined  signal- to-noise  ratio  for  these 
measurements  generally  exceeds  25. 

Purities  of  the  chemicals  used  in  these  experiments  are  as  follows: 
sodium  metal,  99.95%  (Alfa);  helium,  99.995%  (Northeast  Cryogenics);  hydrogen, 
99.995%  (Air  Products);  nitrogen,  99.998%  (Northeast  Cryogenics);  and  hydrogen 
chloride,  99.99%  (Northeast  Cryogenics).  Hydrogen  peroxide,  obtained  as  a  90% 
(by  weight)  solution  from  FMC  Corporation,  Is  purified  by  vacuum 
distillation.  Titration  with  potassium  permanganate  indicates  a  resulting 
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purity  on  the  order  of  99%.  This  corresponds  to  a  purity  in  the  vapor  phase 
of  “93%,  with  the  balance  being  water  vapor. 

Rate  measurements  are  made  using  one  reactant  in  excess  of  the  other 
(detected)  reactant,  thus  ensuring  pseudo-first  order  kinetic  conditions. 
Reaction  times  are  varied  by  changing  the  injector  position  at  fixed  total 
flow  velocity  and  pressure.  Corrections  for  both  axial  and  radial  diffusion 
and  wall  removal  are  made  using  the  procedure  outlined  by  Brown.  This 
method  is  based  on  a  numerical  solution  of  the  equations  describing  diffusion 
and  reaction  in  a  flow  tube.  It  assumes  that  Polseullle  flow  exists  and 
provides  ranges  for  k,  (wall  removal),  and  D  (diffusion  coefficient)  for 
which  the  solutions  have  been  shown  to  be  valid.  Wall  removal  rates  of  the 
measured  species  are  determined  in  separate  experiments  by  varying  the  oven 
position.  The  observed  wall  removal  rate,  as  well  as  the  reaction  rate.  Is 
corrected  for  diffusion  effects.  Diffusion  constants  for  alkali  atoms  and 
hydroxides  were  obtained  from  wall  removal  measurements  In  the  Instance  where 
the  observed  disappearance  of  the  species  Is  diffusion  limited. At  308  K, 
D(Na-He)  =*  0.53  atm  cm^  s"^,  and  D(NaOH-He)  =  0.47  atm  cm^  s“^. 
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An  important  aspect  of  the  data  analysis  Is  the  ability  to  ensure  that 
secondary  reactions  have  no  effect  on  the  concentration  of  the  species  being 
monitored.  This  could  dramatically  affect  the  accuracy  of  the  rate 
measurements.  The  effects  of  secondary  reactions  are  determined  by  modeling 
the  reactions  occurring  In  the  flow  tube.  This  Is  done  using  the  Aerodyne 
PACKAGE  code,^®  a  kinetic  modeling  program  which  numerically  Integrates  the 
differential  rate  equations  for  a  specified  set  of  reactions.  Backward 
reaction  rates  determined  from  the  JANAF  thermochemical  tables  and  the  forward 
rates  are  Included,  to  ensure  that  accuracy  Is  maintained. 

RESULTS 

Na  H?0?  Reaction 

The  NaOH  formation  rate  was  measured  by  directly  observing  the 
disappearance  of  sodium,  with  H2O2  In  known  excess.  The  reaction  has  two 
exothermic  product  channels, 

Na  +  NaOH  +  OH  AH°  -  -130  ±13  U  mol”^  (6) 


>  NaO  +  H2O  AH°  -  -130  ±42  kJ  mol"^  (6') 

so  what  Is  measured  by  monitoring  Na  disappearance  Is  the  total  reaction 
rate.  However,  one  can  add  excess  CO  (~2  x  10^^  cra"^),  which  rapidly  and 
quantitatively  converts  NaO  to  Na,^^“^® 

NaO  +  CO  -*■  Na(^S)  +  GO2  AH°  »  -259  ±42  U  mol~^  (10) 

-  Na(^P)  +  CO2  AH°  -  -56  ±42  U  mol"^  (10') 

Although  the  value  for  k^g  has  not  been  measured.  Indirect  estimates  from 
Reference  (37)  imply  k^g  >  10“^^  cm^  molecule”^  s~^.  The  large  amount  of  CO 


added  Is  required  to  ensure  that  all  of  the  NaO  cycles  back  to  Na  on  a  time 
scale  much  shorter  than  the  time  required  for  reaction  (6)  to  occur. 

A  rate  measurement  with  CO  present  results  In  the  production  rate  for 
only  the  NaOH  branch.  The  results  of  these  measurements  are  shown  In  Figure 
3.  Experimental  conditions  are  given  in  Table  1.  The  rate  constant  for  both 
channels  Is  kg+gi  *  (6.9  ±3.0)  x  10~^^  cm^  molecule”^  s“^,  and  the  fraction 
in  the  NaOH  product  channel  Is  0.61  ±0.10.  The  major  uncertainty  In  the  rate 
constants  kg  and  kg  •  Is  the  H2O2  concentration,  which  is  not  directly 
measured,  but  Is  obtained  as  indicated  earlier. 

An  Independent  determination  of  the  product  branching  ratio  for  reaction 
(6)  Is  obtained  by  observing  the  hydroxyl  radical  via  laser  Induced 
fluorescence  at  308.6  nm.  When  Na  reacts  with  H2O2,  only  the  NaOH  product 
channel  produces  OH.  However,  when  excess  HCl  Is  added  to  the  NaO  and  NaOH 
product  mixture,  additional  OH,  Is  formed  by  the  reaction 

NaO  +  HCl  ♦  NaCl  +  OH  AH°  -  -134  ±42  kJ  mol"^  (11) 

Thus  a  measurement  of  the  amounts  of  OH  produced  before  and  after  addition  of 
HCl  provides  the  relative  amounts  of  NaOH  and  NaO  originally  formed  from  Na  + 
H2O2.  The  observed  NaOH  product  fraction  by  this  method  Is  0.60  ±0.10,  In 
excellent  agreement  with  the  value  obtained  from  the  rate  measurements. 

NaOH  +  HCl  Reaction 

In  the  presence  of  excess  CO,  the  reaction  of  sodium  with  hydrogen 
peroxide  produces  only  NaOH.  Using  this  reaction  as  a  source  for  NaOH,  a 
series  of  rate  measurements  were  made  for  the  reaction  of  NaOH  with  HCl. 
First-order  decays  were  linear  for  more  than  a  factor  of  10  In  fluorescence 
signal,  with  the  HCl  concentrations  ranging  from  1  to  4  x  lo'^^  cra“^.  The 

results  of  these  experiments  are  shown  In  Figure  4,  and  the  rate  constant  for 

this  reaction  Is  k2  *  (2.8  ±0.9)  x  10”^*^  cm^  molecule”^  s'^.  The  uncertainty 
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expressed  Includes  estimated  precision  errors  at  a  95%  confidence  level,  as 
well  as  estimated  errors  In  accuracy. 

Separate  measurements  were  also  made  without  adding  CO.  In  this  case 
NaO,  In  addition  to  NaOH,  Is  present  and  reacts  with  HCl  (reaction  11), 
complicating  the  analysis.  Furthermore,  the  detection  scheme,  based  on 
conversion  to  Na  by  addition  of  atomic  hydrogen,  would  not  appear  to 
distinguish  between  NaO  and  NaOH,  since  Na  atoms  are  also  produced  In  the 
reac  tlon 

NaO  +  H  ->■  Na  +  OH  AH°  »  -155  ±42  U  mol“^  (12) 

r 

Thus,  without  CO,  the  decrease  In  signal  upon  addition  of  HCl  provides  the 
average  of  reactions  (2)  and  (11),  weighted  by  the  branching  ratio  from 
reaction  (6,6').  The  results,  however,  show  the  same  decay  with  and  without 
added  00  (see  Figure  4),  Implying  that  both  NaO  and  NaOH  react  with  HCl  at 
approximately  the  same  rate.  This  observation  Is  reasonable  If  one  considers 
the  hydroxyl  group  on  NaOH  to  act  as  a  quasl-atom  In  respect  to  its  chemical 
behavior  given  Identical  exo thermic 1  ties  for  both  reactions. 

Separate  OH  measurements  confirm  that  atomic  hydrogen  reacts  with  NaO  as 
well  as  with  pure  NaOH.  Although  we  could  not  measure  the  NaO  +  HCl  rate 

constant  directly,  the  observation  that  1)  OH  Is  formed  In  the  NaO  +  H 

reaction  (proving  the  existence  of  NaO),  2)  OH  la  also  formed  upon  addition  of 
HCl  to  NaO,  and  3)  the  decay  rates  with  added  HCl  are  Identical  both  with  or 

without  CO,  Imply  that  NaO  reacts  with  HCl  at  approximately  the  same  rate  as 

NaOH. 


Detector  Corrections  and  Model Ins 


Although  the  plots  of  In  (signal)  vs.  reaction  time  are  linear  over  the 
first  order  of  magnitude  decrease  In  signal,  at  longer  times  they  flatten  out 
at  a  value  typically  a  few  percent  of  the  initial  (zero  reaction  time)  signal 
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(Figure  5).  This  effect  can  be  attributed  to  an  additional  component  to  the 
Na  signal  from  t2ie  reaction 

NaCl  +  H  -*>  Na  +  HCl  -  -21  ±2  kJ  mol"^  (13) 

For  the  case  with  CO  present,  sodium  formed  In  the  detection  region  by 
addition  of  atomic  hydrogen  has  two  sources,  NaOH  and  NaCl,  the  amount  of 
each  depending  on  the  extent  to  which  the  NaOH  +  HCl  +  NaCl  +  H2O  reaction 
has  gone  to  completion.  The  rate  equation  for  the  formation  of  sodium  in  the 
detector  Is 


d[Na] 

dt 


kg  [NaOH]^  [Hi  +  [NaCl]^  [H] 


(14) 


where  the  d  subscript  Indicates  the  concentration  at  the  entrance  to  the 
detector  zone.  Integrating  over  the  detector  reaction  ti.ae  (tj),  assuming 
excess  hydrogen,  we  obtain 


-k  [Hit  -k  [Hit 

[Nal  =>  [NaOHl  .  (1  -  e  ^  ‘^)  +  [NaCll  .  (1  -  e  (15) 

d  d 

If  no  HCl  Is  added,  the  observed  LIF  slgilal  arises  solely  from  NaOH.  We 
have  observed  that  this  signal  (after  correcting  differences  In  diffusion  and 
wall  removal  between  Na  and  NaOH)  Is  identical  to  that  arising  from  only 
sodium  (l.e.,  without  H2O2  added).  This  means  that  within  the  available 
detector  reaction  time,  all  of  the  NaOH  Is  converted  to  sodium,  with  an 
estimated  error  of  20%.  This  sets  a  lower  limit  to  kg  [Hjt^j.  Given  tj  = 

2  ms  and  [H]  ®  2  x  10^^  cm^,  this  requires  that  kg  4  x  10"^^  cm^  molecule"^ 
3“^  ,  The  accuracy  of  t<i  Is  ±0.5  ms.  If  less  than  10%  of  the  Hj  Is 
dissociated  (It  Is  unlikely  to  be  higher  under  our  operating  conditions),  then 
kg  Is  faster  than  the  stated  limit. 
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This  value  for  kg  is  In  reasonable  agreement  with  a  value  of 
I  X  10“^^  10“^^  cm^  molecule'^  s~^  at  300  K)  which  is  used  by 

Hvnes,  et  al.^^  in  fitting  their  flame  data  over  a  wide  range  of  conditions, 
and  also  with  2  x  10"^^  g-800/T  (1,4  ^  10”^^  at  300  K) ,  which  is  used  in 
various  atmospheric  models.  The  former  value  is  claimed  to  be  accurate 

to  within  a  factor  of  two  at  flame  temperatures,  and  the  latter  is  only  an 
estimate.  We  are  planning  to  perform  direct  measurements  of  this  reaction  in 
the  near  future. 

Now  if  we  measure  the  reaction  rate  of  NaOH  +  HCl,  at  short  reaction 
times  the  decay  in  signal  is  logarithmic  because  [NaOHj^j  >  [NaCl]^j,  and 
only  the  first  term  in  Eq.  (15)  is  significant.  At  long  reaction  times, 
[NaCl](j  >  [NaOHl^j  and  the  second  term  becomes  dominant  (Figure  5).  If 
^13  =  kg ,  we  would  never  observe  a  decay  in  [Na]  since  both  terras  in  Eq.  (15) 
would  always  have  the  same  sura.  However,  the  fact  that  an  initial  decay  is 
observed  with  a  later  leveling  of  signal  shows  that  k^g  <<  kg  and,  from  the 
relative  value  of  the  signal  where  it  levels  off,  we  can  estimate  k23. 

To  determine  k^g,  we  modeled  this  reaction  system  using  the  conditions  in 
Table  2  and  rate  set  in  Table  3.  This  rate  set  is  more  extensive  than  the 
above  discussion  implies  because  there  could  be  secondary  reactions  between 
the  excess  H2O2  and  the  sodium  formed  in  the  detector,  as  well  as  other  minor 
OH,  H,  etc.  reactions.  This  model  also  includes  the  undissociated  hydrogen, 
and  the  water  impurity  (~7%)  in  the  H2O2.  We  systematically  varied  kg,  k^g, 
[H],  [H2O2],  and  t^j.  The  results  show  that  all  secondary  reactions  have 
little  effect  (<  5Z)  on  the  calculated  sodium  densities,  and  that  the 
variations  of  Na  with  IH2O2]  and  [H2I  are  small,  in  reasonable  agreement  with 
additional  experiments.  The  value  for  kjg  which  best  fits  the  experimental 
data  is  (8  ±3)/[H],  which  leads  to  value  of  k^g  in  the  range  1  x  to  5  x 

cm^  molecule”^  s~^,  with  a  best  estimate  of  5  x  cm^  molecule”^ 

s“^ .  In  light  of  the  small  exo thermicity  for  this  reaction,  this  value  is  not 
unreasonable.  If  the  rate  constant  is  expressed  in  the  Arrhenius  form, 
k  »  Ae"^/RT^  this  would  correspond  to  only  a  21  kJ  mol"^  barrier,  even  with 
a  gas  kinetic  pre-exponential  term. 
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DISCUSSION 


The  reaction  rate  constant  of  HCl  with  NaOH  Is  found  to  be  in  Its  gas 
kinetic  limit.  Given  the  large  exo thermic Ity,  absence  of  obvious  steric 
effects,  and  the  fact  that  reactant  and  product  states  correlate  on  a  singlet 
potential  energy  surface,  one  might  not  consider  this  surprising.  On  the 
other  hand,  the  strong  covalent  HCl  bond  must  be  broken  during  the  reaction 
with  an  efficiency  approaching  unity.  One  can  calculate  a  rate  constant  for 
this  reaction  using  simple  close  collision  theory.**^  This  procedure  assumes 
that  the  attractive  portion  of  the  Intermolecular  potential  varies  as  Cr“®, 
where  C  depends  on  the  dipole  moments  and  static  polarizabilities  of  the 
reactants,****  and  r  Is  the  Intermolecular  distance.  The  total  effective 
potential  energy  curve  at  large  r  Is  this  attractive  portion  plus  a 
centrifugal  term  which  accounts  for  orbital  angular  momentum.  The  rate  of 
close  collisions  Is  Identified  as  applying  to  collisions  with  total  energy 
sufficient  to  surmount  the  orbital  angular  momentum  barrier  in  this  effective 
potential  energy  curve.  If  It  is  assumed  that  all  close  collisions  react,  and 
that  there  Is  no  activation  energy,  the  rate  constant  calculated  for  NaOH  + 

HCl  Is  1.8  X  10"^®  cm^  molecule”^  ,  In  good  agreement  with  the  measured 
value. 

Few  gas  phase  reactions  of  alkali  molecules  have  been  studied,  none 
Involving  alkali  hydroxides.  The  alkali-hydroxide  bond  Is  Ionic  In  character, 
so  that  (OH)"  Is  a  closed  shell  species,  not  dissimilar  to  Its  Isoelectronic 
and  Isobarlc  analog,  F".  Thus,  In  absence  of  data  on  NaOH,  we  can  use 
available  Information  on  NaF  and  other  alkali  halides  to  try  to  understand  the 
properties  of  the  alkali  hydroxides.  There  have  been  a  number  of  molecular 
beam  scattering  measurements  of  alkali  halides  with  other  alkali  halides  or 
hydrogen  halide  molecules.**^  As  a  class,  they  all  imply  the  existence  of  a 
strong  collision  complex.  In  the  case  of  CsCl  +  KI  -*•  Csl  +  KCl,  the  reaction 
proceeds  without  any  activation  energy  (despite  being  four-centered)  and  has  a 
large  total  cross  section.**®  This  has  been  explained  by  the  large  ionic 
character  of  both  species,  and  by  the  fact  that  alkali  halide  dimers  have 
large  binding  energies  in  the  geometry  of  a  cyclic  planar  rhomboid.**^  This 
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ion-pair  intermediate  formulation  has  also  been  Invoked  in  describing  the  Bal2 
+  HCl  reaction,**^  where  for  large  alkaline  earths  such  as  barium,  we  can  treat 
Bal2  as  (Bal)"^!".  What  we  can  draw  from  these  comparative  systems  is  that  a 
gas  kinetic  rate  constant,  without  any  activation  barrier  is  a  reasonable 
expectation  for  this  system.  In  the  extrapolation  of  the  measured  rate  data 
to  lower  temperatures  characteristic  of  the  stratosphere,  one  would  not  expect 
any  great  dependence  on  temperature  for  the  reaction  rate. 

The  Importance  of  reaction  (2)  in  stratospheric  ozone  chemistry  may  be 
estimated  by  comparing  the  Cl  regeneration  rates  from  HCl  via  reaction  with 
OH,  versus  reaction  with  NaOH  followed  by  photolysis  of  NaCl.  Estimates  of 
stratospheric  NaCl  photolysis  rates  by  Rowland  and  Rogers^**  are  in  the  range 
from  10~  to  10“  s“  .  At  40  km,  the  meteoric  sodium  is  partitioned  among  the 
species  NaOH,  Na02  and  NaO.  Estimates  of  total  stratospheric  sodium 
concentrations  by  Liu  and  Reed^*^  are  on  the  order  of  5  x  10^  cm“^.  if  all  of 
the  meteoric  sodium  were  In  the  form  of  NaOH  at  40  km,  our  measurement  of  k2 
would  give  a  first  order  rate  constant  for  Cl  formation  of  1.4  x  lO"**  s"^, 
providing  that  J^aCl  range  estimated  by  Rowland  and  Rogers.  This 

is  a  factor  of  20  faster  than  the  first  order  rate  constant  for  Cl 
regeneration  from  OH  +  HCl  of  6  x  10”^  s”^,  based  on  an  OH  concentration  of 
10^  cm“^  and  the  value  k5(250  K)  =  5.7  x  10“^^  cm^  molecule”^  Although 

the  partitioning  of  total  sodium  among  NaOH,  NaO,  and  NaO 2  requires  a  more 
detailed  knowledge  of  rate  constants  for  the  processes  shown  in  Figure  1,  this 
calculation  demonstrates  that  even  if  only  5%  is  in  the  form  of  NaOH, 
regeneration  of  Cl  from  HCl  via  alkali  chemistry  would  be  comparable  to 
regeneration  by  the  OH  reaction. 

This  simple  computation  leads  us  to  two  conclusions.  First,  our  results 
for  this  and  the  other  alkali  reactions  described  show  that  meteoric  metal 
reactions  may  have  a  potentially  significant  Impact  on  our  understanding  of 
chemistry  In  the  mesosphere  and  upper  stratosphere.  The  large  values  of  the 
sodium  rate  constants  measured  In  this  study  emphasize  this  possibility. 
Secondly,  It  clearly  motivates  the  need  for  further  Investigation  of  this 
chemistry,  via  detailed  models,  additional  laboratory  measurements  of  reaction 
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rate  constants  and  photolysis  cross  sections,  and  direct  measurement  of 
molecular  sodium  concentration  profiles  In  the  upper  atmosphere. 
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FIGURE  CAPTIONS 


1.  Schematic  diagram  of  atmospheric  sodium  chemistry 

2.  Schematic  view  of  the  flow  tube 

3.  Pseudo-first  order  decays  for  Na  +  H2O2  reaction.  Solid  line  is  total 

reaction  rate;  dashed  line  is  rate  for  branching  to  NaOH  only  (CO  added). 

4.  First  order  reaction  cates  vs.  [HCl]  for  the  reactions  of  NaOH  (NaO)  + 
HCl.  (o)  without  00  added;  (7)  with  CO  added. 

5.  Typical  decay  of  fluorescence  signal  vs.  reaction  time  for  NaOH  +  HCl, 

showing  contribution  to  total  detected  sodium  ( _ )  from  NaOH  ( - )  and 

NaCl  (••••).  Circles  are  experimental  data. 


Table  1  -  Typical  Experimental  Conditions 


Tempera  ture ,  K  ^ 
Flow  Velocity,  m  s“ 
Pressure, , torr 
[Na] ,  cm“  (initial) 
[H2O2],  cm-^ 

[HCl],  cm-" 

[H] ,  cm” 


308 

10.0 


.10 


2.0 

<  10' 

1.3  X  10 
1  -  « 

2  X  10 


13 

10 


12 


Table  2  -  Initial  Conditions  Used  in  Chemical  Model 


Table  3  -  Reaction  Rate  Set  Used  for  Modeling  Na  Chemistry® 


Reac  tlon 

k  (cm^  molecule”^  s~^) 

Reference 

Na  +  H2O2  NaOH  +  OH 

6.9  (-11) 

This  Work 

NaOH  +  HCl  NaCl  +  H2O 

2.8  (-10) 

This  Work 

NaOH  +  H  ->•  Na  +  H2O 

1.0  (-10)  -  1.0  (-14) 

This  Work 

NaCl  +  H  -*•  Na  +  HCl 

1.0  (-10)  -  1.0  (-14) 

This  Work 

Na  +  HO2  -»•  NaOH  +  0 

1.0  (-10) 

Es  tima  ted 

2  HO2  *  H2O2  +  O2 

2.5  (-12) 

39 

0  +  H2O2  -*•  OH  +  HO2 

2.7  (-12)  e'2100/T 

39 

OH  +  H2O2  -*•  H2O  +  HO2 

2.96  (-12) 

39 

H  +  H2O2  -*•  H2  +  HO2 

2.13  (-12) 

40 

H  +  H2O2  -*■  OH  +  H2O 

2.76  (-12) 

40 

H2  +  OH  -►  H  +  H2O 

1.8  (-11)  e~2330/T 

39 

H  +  O2  -*>  0  +  OH 

1 _ 1 q  (—7)  lO.Bbl  p-8316/T 

41 

H2  +  0  -*■  H  +  OH 

3.0  (-14) 

42 

OH  +  HO2  -*•  H2O  +  O2 

4.0  (-11) 

40 

H  +  HO2  ♦  20H 

3.2  (-11) 

40 

H  +  HO2  H2  +  O2 

1.4  (-11) 

40 

H  +  HO2  ♦  H2O  +  0 

9.4  (-13) 

40 

0  +  HO2  +  O2  +  OH 

8.0  (-11)  e~500/T 

40 

20H  -*•  0  +  H2O 

1.0  (-11)  e'550/T 

40 

HCl  +  OH  H2O  +  Cl 

6.6  (-13) 

39 

a)  No  3- body  rates  were  used  because  the  system  Is  at  low  pressure. 
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